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Abstract. The entropic discriminant is a non- negative polynomial associated to a matrix. It 
arises in contexts ranging from statistics and linear programming to singularity theory and algebraic 
geometry. It describes the complex branch locus of the polar map of a real hyperplane arrangement, 
and it vanishes when the equations defining the analytic center of a linear program have a complex 
double root. We study the geometry of the entropic discriminant, and we express its degree in 
terms of the characteristic polynomial of the underlying matroid. Singularities of reciprocal linear 
spaces play a key role. In the corank-one case, the entropic discriminant admits a sum of squares 
representation derived from the discriminant of a characteristic polynomial of a symmetric matrix. 



1. Introduction 

Entropy maximization for log-linear models in statistics leads to the optimization problem 

(1) maximize \x1X2 ■ ■ ■ Xn\ subject to = b. 

Here A is a fixed real d x n-matrix of rank d none of whose columns are zero. The right hand side 
vector b € M'^ is a parameter that is allowed to vary. The problem ([1]) has a unique local solution in 
the interior of each bounded region of the hyperplane arrangement {xi = 0}jg[„] inside the (n — d)- 
dimensional affine space {x G : Ax = b}. The bounded regions are (n — (i)-dimensional convex 
polytopes. Their number, fi{A), is a quantity known in matroid theory as the Mobius invariant. 
The local optima of ([T]) are the analytic centers of these ii{A) polytopes. They are characterized by 

(2) ^4 • X = b and ( — , — , • • • , — ) lies in the row space of A. 

\Xl X2 XnJ 

This translates into a system of polynomial equations in the variables xi, . . . , Xn- It is known [21025) 
that all complex solutions of this system actually lie in M". Thus ij-{A) is the algebraic degree of ([2]). 

The aim of this article is to address the following question: Under what condition on the right 
hand side b do two of the n{A) solutions of polynomial equations represented by ([2]) come together? 
The set of all complex right hand side vectors b G for which this happens is an algebraic variety 
Ha in C^, called the entropic discriminant. Under mild hypotheses on the matrix A, the entropic 
discriminant is a hypersurface and we identify it with its defining polynomial, denoted ff^(b). 
This is a non- negative polynomial whose real zeros lie in certain linear subspaces of codimension 2. 

Example 1. Let d = 3 and n = 5. The following 3x5-matrix has Mobius invariant fJ-{A) = 4: 

A = 




The entropic discriminant of A is a homogeneous polynomial ^/'^(^i, ^2, ^3) of degree 8. It equals 
288blbl{blbl + blbl + blsl + blsl + blsl + blsl + blsl + 6§s§) + 17736|&| + 720&i6i(s?si + 6?s§) 

+ lQ2{blbisl + b\slsl + blblsl + blslsl) + 1216(626252^2 _,. 52^2^2^2 ^ fe2^2^2g2 ^ 52^2^2^2) ^ 2^Qblslslsl 
+320(6262,2,2 ^ ^2^2,2,2 ^ ^2^2,2,2 ^ ^2^2,2,2 ^ ^2^2,2,2 ^ ^2^2,2,2 ^ ^2,2,2,2 ^ 62,2,2,2)^ 

where si = 61 — 62; ^2 = 61 — ^3) and S3 = 61 — 62 — 63- Thus HAih) is a sum of squares of quartics 
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It coincides with the discriminant of the fohowing system of equations in three unknowns: 

1/zi + l/(zi + Z2) + 1/(^1 + ^3) = h, 

1/Z2 + 1/(Z1 +Z2) = 62, 

I/Z3 + 1/(Z1 +Z3) = 63. 

These equations are equivalent to ([2]) if we take {zi, -22, ^3) to be coordinates for the row space of A. 
There are four solutions for any b = (61, 62, ^3) G C^. They are distinct if and only if Ha{^) 7^ 0. 
The entropic discriminant -ff^(b) is a non- negative polynomial having precisely four real zeros: 

(3) Vr{Ha) = {(0:1:0), (0:0:1), (1:1:0), (1:0:1)} C P^. 

The complex variety Vc{Ha) is a curve of degree 8 in the projective plane with coordinates (61:62:63). 
That curve is singular at its four real points. In addition, it has 16 isolated complex singularities. o 

We shall study the systems ^ for arbitrary d, n, and A. The following is our main result: 

Theorem 2. Let A be a real d x n-matrix of rank d whose columns span > d + 1 distinct lines. 
The entropic discriminant is a hypersurface, defined by a homogeneous polynomial HaO^) of degree 

(4) degHAih) = 2i-lf-idxi0)+x'm, 

where x{t) is the characteristic polynomial of the rank d matroid of A. For generic matrices A, 
this degree equals 2(n — d){^^. The polynomial -ff^(b) is non-negative for all arguments in W^. 

The generic degree 2(n — d)(]J~2) is always an upper bound on the degree of the entropic dis- 
criminant, and equality holds when the matroid of A is uniform; cf. Proposition [33l For example, 
for generic matrices A of size 3x5, the degree of i?A(b) equals 16, and not 8 as in Example [TJ 

This article is organized as follows. In Section [2] we examine the polar map of a product of linear 
forms. The entropic discriminant is shown to coincide with the branch locus of that polar map. 
For example, consider the polar map of the binary form /(^i, -22) = •2i(-Zi + 2z2){zi + 'iz2){zi +0^2): 

: ^ Pi, (zi :Z2) ^ (|^(^i'^2) : ^(^1,^2))- 

The branch locus of this map consists of the four zeros of the binary quartic HA{bi, 62) in Example [3] 
below. This connects our study of HA{h) to the topological theory of hyperplane arrangements [HE], 
and to topics in classical algebraic geometry that are found in Chapter 1 of Dolgachev's book [5]. 

Section [3] is concerned with the important special case n = d+1. Here the entropic discriminant 
has expected degree d{d — 1) and we can write it explicitly as a sum of squares. This expression 
is derived from known results on the discriminant of the characteristic polynomial of a symmetric 
matrix [2l ll3^ [n]. We then apply this to resolve two problems left open in the literature, namely 
Conjecture 1 in [T] concerning the discriminant of the derivative of a univariate polynomial, and 
Conjecture 7.9 in [23j concerning the real critical double eigenvalues of a net of symmetric matrices. 

For any linear subspace C of C", its reciprocal is defined as the Zariski closure of the set 

(5) {(irr'ii;'---'i^)^'C" : (^i,^2,...,n„)G/:n(CT}. 

In Section [5] we study the geometry of the reciprocal plane C~^, further extending the line of the 
work from Proudfoot-Speyer [19] to Huh-Katz [12j. We identify a minimal system of defining 
equations for C~^, we characterize the singular locus of and we determine all tangent cones. 
The relationship between that singular locus, the ramification locus of the map A : C^^ — )• P'^^^, 
and the entropic discriminant HA{i>) is studied in detail in Section [71 In Corollary [37] we show that 
the real variety defined by the polynomial //^(b) is a union of linear spaces of codimension 2 in 
p*^-!. We saw this already for one instance in Example [U where d = 3 and the real variety is finite. 
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Theorem [2] is proved in Section [6l However, one subtle but essential point needs to be taken 
care of before that proof. In order for Q to be the correct degree, a more refined notion of 
entropic discriminant is required. Namely, we shall define Hyi(h) as the polynomial defining the 
cycle-theoretic branch locus of the restriction to C^^ of the linear map A : C" — > C^, where £, is 
the row space of A. The following example justifies this "fine print" in Definition | 



Example 3. Let d = 2, n = 4 and A = ^ ^ ^| where a is a real parameter. For general 
values of a, the entropic discriminant is irreducible and has degree 4, as predicted by Theorem [2j 

HA{bi,b2) = (22680-* - 97200^ + 11664a2)fe4 _ (aoOOa* - 12528a3 + 12960a2 + 5184a)6f62 
+ (1744a4 - 7980a3 + 10584a2 - 2160a + 5184)6261 



-(SOOa^ - 2612a3 + 46800^ - 3888a + 4320)6i6^ + (eSa^ - lOOa^ + 999a2 - 1350a + 1188)6: 



For special values of the parameter a, this expression factors over Q. For a = 6, it is the square 
972(366^ - 246162 + 56|)2. Thus, here the four points of Vc{Ha) in are two double points. o 

Our initial motivation for embarking on this project was a model in theoretical neuroscience 
proposed by Hillar and Wibisono [lOj. These authors investigate the retina equations which char- 
acterize the maximum entropy distribution for a graphical model G with n edges having continuous 
random variables on d nodes that represent the firing pattern of d neurons. Their equations are 

(6) Yl " for i = l,2,...,d, 

where Af{i) is the set of all nodes that are adjacent to the node i. The real numbers 61, 62, . . . , 6^ 
are parameters that serve as the sufficient statistics of the desired maximum entropy distribution. 

To fit the system ([6|) into our framework, we introduce new unknowns Xij = l/{zi + Zj) for all 
edges {i,j} € E[G). This translates ^ into the linear system ^-x = b, where A is the node-edge 
incidence matrix of G and x = {xij : {i,j} € E{G)) is a column vector of unknowns. Of course, 
these unknowns obey the additional constraints that x must lie in the reciprocal plane C^^, where 
C is the row space of A. Thus the retina equations of Hillar and Wibisono fit our format ([2]): 

(7) ^-x = b and X E C'K 

The entropic discriminant Hji^{h) characterizes measurements b for which the retina equations ^ 
or ([7]) have multiple roots. Of particular interest is the case n = (2) , when G = is the complete 
graph, and the sum in ^ is over j € {1, . . . ,n}\{i}. The characteristic polynomial Xd{t) of the 
corresponding matroid was computed by Zaslavsky |27j . in his work of colorings of signed graphs: 



(8) xs) = E({'} + ^{V})(*-i 



(2) 
k ■ 

k=0 



Here {^} is the Stirling number of the second kind and (x)^^^-^ = x{x — 2) ■ ■ ■ {x — 2k) is the 
generalized falling factorial. One can also compute Xdit) with the exponential generating function 

(9) T.Xd{t)-^_ = (l + x)•(2•exp(x)-l)(*-^)/^ 

found in [22\ Exercise 5.25]. Using these formulas, one obtains the first few values of the degree of 
i?A(b) and of the number of solutions of the retina equations on the complete graph G = Kd'. 

d = 4567 8 9 10 

(10) deg(i7A(b)) = 22 270 3148 38990 524858 7705572 123087958 

H{A) = 7 51 431 4208 46824 586141 8161237 
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The requisite combinatorics is developed in Section HI It covers material from matroid theory, 
focusing on geometric interpretations of the characteristic polynomial and the Mobius invariant. 

2. The polar map of a product of linear forms 
The d X n-matrix A = (aij) determines a product of linear forms in d unknowns z = (zi, . . . , Zd)' 

n d 

(11) /(Z) = \{{Y.^^3Zi). 

j=l i=l 

The hypersurface Vc{f) is an arrangement of n hyperplanes in the complex projective space P'^"^. 
The polar map of this hypersurface is the rational map 

The base locus of V/ is the singular locus of Vc(/), and this is the union of all codimension-2 
strata in the hyperplane arrangement. If the columns of A are linearly independent then V/ is 
the Cremona transformation of classical algebraic geometry, and, in general, the polar map V/ is 
also known as the polar Cremona transformation [6]. The Jacobian of Vj is the Hessian of the 
polynomial /, that is, the symmetric matrix of second derivatives. We consider its determinant 

Hess(/) = det ' ^ 



This is a homogeneous polynomial of degree d{n — 2). Its zero set in F'^-^, denoted by yc(Hess(/)), 
is also referred to as the Hessian of /. We are interested in the image of that hypersurface under V/. 

Proposition 4. The entropic discriminant equals the image of the Hessian under the polar map: 

(12) Vc{Ha) = closure of Vf{Vc{Ress{f))\Vc{f)). 

Proof. Let denote the reciprocal of the subspace C spanned by the rows of A, regarded as a 
subvariety of P"~^. The variety is the closure of the image of the map P"'^^ --4 that takes 
a general point z = (zi : • • • : z^) in P^^^ to (zA)"^ = ((XliLi Oji^j)""^ : ' ' ' : iYli=i am^:^)"^) in 
pn-i^ The polar map is the composition of this map with the linear projection p*^"! --4 P"'"^, x 1— > 
Ax. In symbols, we have V/(z) = A(^{zA)~^^ . This observation shows that the fiber of Vj over a 
general real point b G Im(Vj) consists of IJ-{A) real points in P"'"^, namely, the points represented 
by the analytic centers in the arrangement defined by the coordinate hyperplanes in the affine space 
{x G M" : Ax = b}. This result was also obtained by Dimca and Papadima in [H Corollary 4 (1)]. 

For special complex points b G P*^"^, two of its fJ-{A) preimages under V j may coincide. At such 
a preimage z of multiplicity > 2, the Jacobian of drops rank, and the Hessian of / vanishes at z. 
Conversely, points z outside the hyperplane arrangement Vc{f) at which the polynomial Hess(/) 
vanishes must be double roots of the system of equations Vj(z) = b. Since the parametrization 
z I— 7- x = {zA)~^ maps P'^^^ birationally onto the reciprocal plane C~^, such double roots appear if 
and only the intersection fl {x G P"~^ : Ax. = b, X1X2 • • • x„ 7^ 0} has a point of multiplicity 
> 2. This condition on b is the geometric definition of the entropic discriminant Hj^. □ 

We have not yet addressed the question whether the entropic discriminant actually has codi- 
mension 1, and this may in fact not be the case. For instance, if A is the identity matrix 
and / = Z1Z2 ■ ■ ■ -Zd then is the classical Cremona transformation on P<^~1 and Hess(/) = 
(— 1)/*^"^. Here, the Hessian coincides with the hyperplane arrangement, and the entropic 
discriminant is not a hypersurface. We shall see that this is essentially the only exceptional case. 
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The matrix A = (aij) is called basic if its column rays lie on d distinct lines in W^. Since A has 
rank d and no zero columns, this means that the distinct column directions form a basis of M*^. 

Corollary 5. If A is not basic then the entropic discriminant is a hypersurface in P'^^^. 
Proof. A classical formula [15^ p. 660, Ex. 10] for the Hessian determinant of / states that 

(13) Hess(/) = (n-l) f ''-^ ■ det{Aif H {aikZi + a2kZ2 + --- + adkZd? 

where Aj denotes the dxd-submatrix of A with column indices /. If A is not basic then there are 
at least two linearly independent summands. This implies that the Hessian hypersurface is not 
contained in the hyperplane arrangement Vc(/). The polar map is a finite-to-one morphism on 
the open set P^-^\Vc(/), and hence it maps the Hessian to a hypersurface in P"^ namely -ff a- D 

Corollary 6. For any non-basic A, the polynomial //^(b) is homogeneous and nonnegative on M*^. 

Proof. It suffices to prove this for the square-free polynomial iJyi(b) that vanishes on Vc{Ha)- 
Indeed, if HA_{h) is homogeneous and nonnegative then so is any real product of its factors. 

Homogeneity is straightforward since the geometric definition ensures that b € Vc{Ha_) implies 
Ab € Vc{Ha)- To show non- negativity, let K denote the subfield of M which is generated by the 
entries of A. We regard the entries of b = (bi, . . . ,bd) as indeterminates over K. Let L be the 
algebraic closure of the rational function field K(bi, . . . ,bd). Then the equation V/(z) = b has 
fJ.{A) distinct solutions with coordinates in L. We substitute these solutions into the sum in (jl3p 
and we take their product in the field L. The result is a sum of squares in L that is a symmetric 
polynomial in the roots. It is invariant under the action of the Galois group of L over K(bi, . . . , bd) 
and thus lies in K{bi, . . . , bd). The sum of squares representation over L ensures that this rational 
function is non-negative under all specializations of b to M at which it does not have a pole. The 
numerator of this rational function is a product of the factors of i?y4(b). We conclude that HaO^) 
does not change signs on M.'^. Hence, either iJyi(b) or —HaO^) is non-negative on W^. □ 

The above argument shows that Hyiih) is non- negative but it does not furnish a representation 
of HA{h) as a sum of squares of polynomials. We also note that the computation of HA{i>) from 
Hess(/) is a task of elimination theory that is quite non-trivial even for moderate values of d and n. 

One case where the elimination problem can be solved more easily is d = 2. Here f{zi,Z2) is 
a binary form of degree n enjoying the property that all its zeros on the line are defined over 
M. The polar map V/ takes the complex projective line P^ to itself. This map has degree n — 1, 
i.e. the fiber over a general point b € P^ consists of n — 1 points. We are interested in those points 
b on the line P^ for which two or more of the points in its fiber collide. The Hessian of / equals 

[l-n) ■ ^ {aua2j - aija2if ■ (aifc^i + a2kZ2) 

l<i<j<n keln]\{i,j} 

This is a binary form of degree 2n — 4, so it defines a configuration of 2n — 4 points in P^. All points 
have non-real coordinates. The entropic discriminant of / is the image of these 2n — 4 points under 
the polar map V/. Proposition H] gives the following rule for computing the entropic discriminant: 

df df 

(14) HAibiM) = Resultantz(Hess(/(z)),62^(z)-6i^(z)). 

This formula can be rewritten as the discriminant of a binary form: 

df df 

(15) HA(bi,b2) = Discriminantz(&2T^(z) - ^17^ — (z)). 

OZi OZ2 




\2 
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The binary form HAibi, 62) has degree 2n — 4 provided no two columns of A are parallel. Being 
nonnegative, the entropic discriminant is a sum of squares of binary forms of degree n — 2 over M. 

Example 7. Let n = 3 and consider a general binary cubic with real zeros: 

/ = {aiiZi + a2iZ2)iai2Zi + a22Z2){ai3Zi + a23Z2)- 
The sum of squares representation in (fT3]) tells us that the Hessian of / equals 

{aiia22-a2iai2f{ai3Zi+a23Z2f+iaua23-a2iai3)'^{ai2Zi+a22Z2f+iai2a23-a22ai3)'^{auz 

For any invertible matrix U, the entropic discriminant satisfies HjjA{Uh) = ff^(b). This implies 
that /7^(b) can be written in the 2x2 minors pij of the matrix 2 x 4-matrix {A,h). We have 

(16) HA(h) = {pi2-P34:f + {Pl2.-P2if + (P23-Pl4)^- 

For n = 4, an expression for i7/i(b) in terms of the 2x2 minors the 2 x 5-matrix {A, b) is 

(--^^X b?2P34P35P45)^ + (P?3P24P25P45)^ + (PuP23P25P35)^ + (^14^13^15^45)^ + {VlZPlAPl^P'i^f 

+ {Pl2PliPl5P25f + i(P23P24P34P?5)^ + liPliPUPiiPhf + | (Pl2Pl4P24Pi5)^ + 5(^12^13^23^45)^- 

At present we do not know how to extend the formulas (|16p and (|17p to n > 5. o 



It is natural to ask how the formulas (|14p and (jl5p would generalize to d > 3, and the answer is 
given by the projective duality between the entropic discriminant and the Steinerian hypersurface [U 
§1.1.6]. If / is any homogeneous polynomial of degree n in z = (zi, . . . , z^) then its Steinerian is 

/ df df df \ 

(18) St/(ci,C2,...,Crf) = Discriminantzl ci — (z) + C2 — (z) H hQ— (z)j. 

In this formula, we are taking the discriminant of a form of degree n — 1, namely, the polar of / 
with respect to a generic point c. Corollary 1.2.2 in [5] tells us that the hypersurface defined by 
Stj(c) is dual to the image of the hypersurface defined by IIess(/(z)) under the polar map Vj. 

In our situation, the given form / is a product of linear forms as in (lllh . and some care needs to 
be taken in removing contributions from singularities. Indeed, the Steinerian Stj of a hyperplane 
arrangement is supported on that same hyperplane arrangement plus an extra component. It is 
this extra component we are interested in. We call this hypersurface the residual Steinerian of /. 

Corollary 8. The entropic discriminant of a dxn-matrix A is the hypersurface inF'^^^ projectively 
dual to the residual Steinerian of the arrangement of n hyperplanes given by the columns of A. 

Let us briefly revisit the case d = 2 from this point of view. We saw that the entropic discriminant 
consists of 2n — 4 points on a projective line with coordinates (bi : 62)- The Steinerian consists of 
2n — 4 points on the dual projective line with coordinates (ci : C2). In our formulas (I14p and (jlSp 
we tacitly identified these two lines and their point configurations via (ci : C2) = (—62 : 61). 

For d > 3, the formula (jlSp is less useful for the purpose of computing HaO^) because dualizing 
the residual Steinerian in a computer algebra system is hard. Instead, we find it preferable to use 

(19) {HA{h)) = ((Hess(/(z))) + (2x2-minors of the 2xd-matrix(b,V/))) : (V/)°°. 
This ideal-theoretic reformulation of (I12p is the direct generalization of (jl4p to d > 3. 

Nevertheless, the (residual) Steinerian of a hyperplane arrangement remains a beautiful topic in 
geometry, and its interplay with the combinatorics of the entropic discriminant certainly deserves 
further study. We close this section with an illustration of this for lines in the plane P^. 

Example 9. This example was worked out with help from Igor Dolgachev. Let d = 3 and suppose 
the matroid of A is uniform. Thus Vc(/) is an arrangement of n lines in general position in P^. 
By Theorem [21 the entropic discriminant Ha is a curve of degree 2(n — l)(n — 3). Its singular 
locus consists of the n columns of A. By dualizing, we obtain the Steinerian St/, a curve of degree 
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3(n — 2)^. Each of the n hnes occurs with multiplicity n — 2 in the Steinerian. Removing these lines, 
we find that the residual Steinerian is a curve of degree 3(n — 2)^ — n(n — 2) = 2(n — 2) (n — 3) . o 

3. The codimension-1 case 



The discriminant of the characteristic polynomial of a symmetric matrix is non-negative because 
real symmetric matrices have only real eigenvalues. The study of this discriminant is a classical 
subject in mathematics, going back to an 1846 paper by Borchart [2J. Explicit representations of 
this discriminant as a sum of squares were also presented in work of Newell |17j , Ilyushechin [13] , 
and Lax |14j . See [231 §7.5] for an exposition, and work of Domokos [7] for the state of the art. 

In this section we establish a relationship between this subject and the entropic discriminant. 
We focus on the case n = and we express ffyi(b) as a specialization of the discriminant of the 

characteristic polynomial of a symmetric matrix. We shall use this to derive the following result. 

Theorem 10. Let A he a non-basic matrix with d rows and n = d+1 columns. Then the entropic 
discriminant HaO^) is a sum of squares of polynomials. Moreover, if the entries of A are rational 
numbers then -ffA(b) is a sum of squares in Q[bi, . . . ,bd\. 



Example 11. If d = 3, n = 4 and A 




then HA{bi, 62, ^3) is the sum of 10 squares 



lhi{b2 



-&3)^ + i( 



n- 



+Ubib2+bibs-2b2b3)Hbi 
-761 



b2fblbl + ^(56162-96163-1462+186263)^6? + ^(56162-36163-862+66263)26? 



262 



+ 



T_ 

108 



(56162+36163 



-262- 



-6626: 



'203; 



^62 + 



1 
216 



(I36162-2I6163 



126263 + 2762)26? + -^(56? 62 - 76? 63 - 76i62 + 4616263 + 96i6i + Ub^ - 18626^2 



+ 216(5^1^2 - 2I6163 + 62 - 126263 + 2762)26? + 3^(56262 - 6?63 - 46i62 - 8616263 + 86^63^2 
This expression is derived from the sum of 10 squares found at the top of page 97 in [23 



Proof of Theorem\l(A Let ^ be a non-basic d x (fi+l)-matrix and let v S W^^^ span the kernel of 
A. If V has a zero coordinate, say Vd+i = 0, then we can reduce our analysis to a smaller case, 
namely, a {d—1) x d-matrix obtained by taking the columns of A modulo the last column. Hence 
we may assume that all coordinates of v are non-zero. 

Next, we claim that it suffices to prove our assertions for the special case where 



and w = (1,1,1,... ,1)^ 



That this suffices is ensured by the following transformation rule for the entropic discriminant: 







(I 





• • 


• 


-1\ 









1 


• • 


• 


-1 


(20) 


A = 








1 • • 


• 


-1 











•• 


• 1 


-1/ 



(21) 



HuAD{h) = HAiU-'h 



This identity holds for any invertible dxd-matrix U and any invertible diagonal nxn-matrix D, 
and its validity is easily seen from the geometric definition of Ha. We here use this for n = d + 1. 

We now fix A and v as in ([20]) . Then C = rowspace(j4) is the hyperplane xi + X2 + ■ ■ ■ + Xn = 0. 
Its reciprocal is the hypersurface of degree d in P"^ that is defined by the polynomial 



/XI + 



(22) 



En-. 

1=1 



det 



X2 + Xr, 



\ 
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This symmetric determinantal representation of the (n — l)st elementary symmetric polynomial is 
taken from [20j . The linear system = b is equivalent to 

(23) Xi = hi+Xn for i = 1, 2, . . . , n — 1. 

Thus the points satisfying ^ can be computed by substituting (j23p into (j22p and equating the 
resulting univariate polynomial to zero. Setting t = Xn, the solutions to ([2]) correspond to zeros of 









1 


1 • 


■ A 






1 


2 


1 • 


• 1 


(24) 


Phii) = det(ti? + diag(6i, 62, . . . , 6(i) ) where E = 


1 


1 


2 • 


• 1 






V 


1 


1 •• 


• y 



In particular, i?^(b) equals the discriminant of the univariate polynomial Pbit)- The following 
proposition applied to E and X = — diag(6i, . . . ^bd) completes the proof of the theorem. □ 

Proposition 12. Let E G 5g symmetric positive definite matrix and X a symmetric matrix 

of indeterminates. Then the discriminant of the generalized characteristic polynomial det(tE — X) 
with respect to t is a sum of squares in Q{Eij ■ 1 < i, j < m)[Xij ■ 1 < i, j < m]. 

Proof. Since E has a Cholesky factorization E = MM'^, it follows that 
(25) det{tE-X) = det{E) ■ det{tl - M'^XM-^). 

We get a sums of squares formula from the known representations of the discriminant of the 
characteristic polynomial of a real symmetric matrix. However, our emphasis lies on the rationality 
of the desired formula. Following |23t §7.5], let X = M~^XM~'^ and consider the linear map 

^ SymaR™ , Z ^ [X,Z] = XZ - ZX 

that takes a skew-symmetric matrix to the commutator with X . 

Let {Wij = Ci A Cj : 1 < i < j < m} be the standard basis of skew-symmetric matrices and 
likewise {Sij = Ci ■ ej : 1 < i < j < m} the standard basis of symmetric matrices. Let $ be the 
("^2"^) ^ (^) -matrix representing the linear map in the chosen bases. By choosing suitable bases, 
it can be seen that the eigenvalues of are the squared pairwise differences of the eigenvalues 
of X. Hence the determinant of is the discriminant of det(t/ — X). The sum of squares 

representation can be obtained by applying the Binet-Cauchy theorem. 

To get a rational representation we apply the above reasoning to the slightly altered map 

Z ^ [Z,X]e := E-^XZ - ZXE-^. 

It is clear that a representation in the standard basis is over <Q{Eij) and hence yields an appropriate 
sum of squares. To see that this actually yields the discriminant for the generalized characteristic 
polynomial, choose bases W^j = M~'^ and S[j = M~'^ Sij and verify 

[X,W^^]e = M~^[X,Wij]M-\ 
Hence, a representation in the new bases is given by <I? above. □ 

Evaluating the discriminant of p^ (t) in ()24p leads to the following data concerning the monomial 
expansion of the entropic discriminant HA{h) of the particular matrix A in (I20p : 



(26) 



d 


2 


3 


4 


5 


6 


degree of HA{h) 


2 


6 


12 


20 


30 


number of monomials 


3 


19 


201 


3081 


62683 


leading (lex) monomial 


bl 


bfbl 


bm 


bfblbjbl 


bl%lblbibl 
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Remark 13. The entropic discriminant HA_{h) is a symmetric polynomial in 61, ... ,6^ since the 
set of rows of the matrix A is invariant under permutations. It thus admits a unique representation 
as a polynomial in the elementary symmetric polynomials 

Cfc = X] bi^bi.^ ■ ■ ■ bi^ for A; = 1,2,... ,d. 

l<.ii<---<if^<n 

For example, if d = 4 then the 201 terms in 61, &3) ^4 translate into only 16 terms in ei, 62, 63, 64: 

HAih) = 432efe| - 432efe2e3e4 + USefel + 108e?e^e4 - 36e?e^e^ - 2160efe2e| + 1800eie^e3e4 
+ 120e?e^e4 - 540eie2ej^ - 405e^e4 + 135e^e^ + 2400eie3e| + 1800eie| - 2700e2e^e4 + 675e| - 2000e|. 

As an application of our theory, we are now able to answer two questions from the literature. 
The first deals with the discriminant of the derivative of a univariate polynomial. According to 
Alexandersson and Shapiro [T] Conjecture 1], Frank Sottile and Eugene Mukhin formulated this 
conjecture at the AIM meeting "Algebraic systems with only real solutions" in October 2010. 

Corollary 14. The discriminant of the derivative of a univariate polynomial f{t) of degree n is a 
sum of squares of polynomials in the differences of the roots of f{t). 

Proof. Let P„ = discrf (/'(t)) . We shall write P„ specialization of the entropic discriminant 
and use the sum of squares decomposition given in TheoremllOl Consider the univariate polynomial 
f{t) = rir=i(* ~ 0,1). Notice that Xi = t — ai provides a parametrization for the one-dimensional 
affine space {^x = b}, where we take A as in (j20p and bi = a„ — Oj for i = l,...,n — 1. We plug 
this parametrization into the polynomial (j22p that defines This yields the derivative f'{t) = 
Sj=i Wijkj{'t~ '^i)- Thus f'{t) equals the polynomial pb(0 of dMl) whose discriminant (with respect 
to t) equals ///^(b). We conclude that T>n equals the entropic discriminant HaHou — flj)ie[n-i])- 
Using Theorem [TOl we conclude that Vn is a sum of squares in Q[b] = Q[(a„ — Oj) : i G [n — 1]]. □ 

Our techniques can also be applied to answer a question that was left open in [231 §7-5]. Namely, 
we conclude this section by proving Conjecture 7.9 of [23,]. 

Corollary 15. There exist three real symmetric d x d-matrices Co,Ci,C2 such that all (^^3^) pairs 
of complex numbers (x, y) for which Cq + xCi + 7/C2 has a critical double eigenvalue are real. 

Proof. Consider the symmetric matrix X = M^^XM^'^ with X = diag(6i, . . . ,6^) in the proof of 
Proposition [T2j Its entries are linear forms in b = (61,621 • • • ^bd)- We replace the unknowns bi by 
generic real affine-linear forms in two variables x and y, say bi = Wi + UiX + Viy for i = 1, . . . ,d. 
The symmetric matrix resulting from this substitution is a net of real symmetric d x d-matrices: 

X = Co + xCi+yC2. 

The real values of (x, y) for which this matrix has a critical double eigenvalue corresponds to the 
intersections of this affine plane with the real variety of HaO^), with A given in (|20p 

We claim that the real radical of the entropic discriminant is the codimension 2 ideal 
(27) vTMb)) = n ^ n {bi-bj,b,-bk). 

'^<i<j<d l<i<j<k<d 

This identity follows from the geometric description of //^(b) in terms of colliding analytic centers. 
Indeed, the hyperplane arrangement defined by {xi = 0} in {x : Ax = b} consists of n = d + 1 
points on a line. They form d bounded segments. The analytic centers of two segments collide 
if and only if three of the d + 1 points coincide. There are such C^^^) triples, each imposing a 
condition of codimension 2. They are expressed by the ('^^^) prime ideals in the intersection (127p . 
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Since the real variety of HA{h) is a union of C^^^) real linear spaces, each of its intersection 
points with the plane b = w + ux + vy is also defined over the reals. Therefore all C^^^) symmetric 
matrices with a critical double eigenvalue in the net Co + xCi + yC2 have real entries. □ 

Remark 16. The validity of the identity (]27p above rests formally on the inclusion Sing(>C~^) C 
TZa- This inclusion will be proved in Example [32] for the particular matrix A in ()20p . The issue of 
determining ^ {ifyi(b)) in general, for n > d + 2, is very subtle. This will be our topic in Section [71 

4. Matroids and Graphs 

In this section we discuss the notions from matroid theory which are needed for the statement 
and proof of Theorem [2j We also discuss various matroids arising from graphs, including those 
representing the Hiller-Wibisono model ([6|). Matroid theory is a classical subject in combinatorics 
with many (axiomatic) paths leading to it. For us, matroids come in the form of matrices and 
hence we take the concrete approach via realizable matroids. For more on this subject see |18ll22j . 

Our given matrix A = {Ai,A2, . . . , An) G R*^^" is identified with an ordered collection of n vectors 
that span d-space. The corresponding matroid M = M(A) records all linear dependencies among 
these vectors. A subset / C [n] = {1, 2, . . . , n} is called (in) dependent whenever Aj = {Ai : i G I) 
is linearly (in) dependent. The rank rk(/) of / is the rank of Aj. The rank of the matroid M is the 
rank of A. A circuit is an inclusion-minimal dependent subset, and / is independent if it does not 
contain a circuit. A subset F C [n] is a flat if F equals {i G [n] : Ai G span{AF)}, that is, if F 
precisely indexes a collection of vectors contained in some linear subspace. Equivalently, F is a flat 
if and only if it meets every circuit in > 2 elements or not at all. Flats will play an important role 
in Sections [5] and [H The lattice of flats L{M) is the collection of all flats, ordered by inclusion, 
with minimal element = {z : >lj = 0} and maximal element 1 = [n]. The lattice of flats represents 
combinatorial information about the containment relations of the various subspaces spanned by 
subsets of columns of A. It is one of the central objects in the enumerative theory of matroids. 

Associated to L = L[M) is its Mobius function fi^ : Lx L ^ "Z, which is defined by fJ^iiF, F) = 1, 

f^L{F,H) = - Yl t'LiF^G) 

F<G<H 

\i F < H, and ^l{F, H) = otherwise. The characteristic polynomial of M is defined by 

XM{t) = /^L(M)(6,F)t^^'(^-^)-'-'=(^). 
FgL(A/) 

The (unsigned) Mobius invariant of M, or of the matrix A, is the positive integer 

fi{A) = ix{M{A)) = I^L(6,i)I = (-I)'^xm(O). 

Evaluations of the characteristic polynomial have nice combinatorial interpretations in terms of 
hyperplane arrangements [22]. For a generic b G M"^, the set of n hyperplanes {xi = 0}jg[„] induces 
an affine hyperplane arrangement in {x : ylx = b} = M""*^. The Mobius invariant /i(^) equals 
the number of bounded regions of this arrangement. This fact played an important role in [3| §3]. 
The proof is a straightforward deletion-contraction argument, using that ^{A) and the number of 
bounded regions in {Ax. = b} adhere to the same recurrence relations. This number is independent 
of the choice of generic b and is also known as the beta invariant of the extended matrix (^,b). 
The proof of the following observation illustrates the typical line of arguments in matroid theory. 

Proposition 17. The Mobius invariant ^i{A) equals 1 if and only if the matrix A is basic (defined 
in Corollary \^ if and only if its geometric lattice L{M) is the Boolean lattice of all subsets of [d]. 
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Proof. By the definition of the lattice of flats, we can assume that A has no zero columns and that 
no two columns are proportional. Hence, the matrix A is basic if and only if M(A) is isomorphic 
to the uniform matroid U^^d whose lattice of flats is the Boolean lattice of all subsets of [d]. The 
Mobius invariant of the matroid Ud^d is = 1; see Example 1181 below. Conversely, if A is non- 
basic, there is a column e S [n] that is not an isthmus, that is, not contained in every basis. For 
such an element e, the Mobius invariant satisfies the deletion-contraction identity 

fi{A) = fi{A\e) + fi{A/e). 

By Rota's Sign Theorem, the Mobius invariant is always a positive integer and hence ^{A) > 2. □ 

We have now defined the combinatorial ingredients for the degree (HD of the entropic discriminant. 
With this in place, we derive the value of that degree for generic matrices A stated in Theorem [2j 

Example 18 (Uniform matroids). A generic d x n-matrix A with d < n gives rise to the uniform 
matroid M = Ud,n in which every subset of cardinality < d is independent. The corresponding 
lattice of flats is a truncated Boolean lattice in which a subset F C [n] is a flat if and only if < d 
or F = [n]. The Mobius function on the Boolean lattice is fi{F,G) = (— for F <^G. Hence 

Xu,,St) = t'-nt'-' + ••• + + (-1)'^^:^). 

Note that t = 1 is always a zero of the characteristic polynomial. The number of solutions of the 
equations ([2]) for generic A equals ij-{A) = {^Zl)- The degree (|4]) of the entropic discriminant equals 



2i-l)''-idxu,JO) + x'u,„m = 2 



n — 1\ ( n 

d\ 



d-l \d-l 



2(n - d) 



n — 1 
d-2 



As we will see in Proposition 1331 below, this quantity is an upper bound for fixed n and d. o 

Graphical matroids are an important class of examples. Let G be a graph on d nodes with n 
edges and c connected components. For an arbitrary but fixed orientation of the edges, let Aq be 
the d X n incidence matrix of node-edge pairs, with entries — 1, if the node is in-coming, out- 
going, or non-incident for the edge. Reorienting an edge of G results in scaling the corresponding 
column of Aq by —1 and hence leaves the matroid Mq = M{Ag) invariant. Note that Aq has 
rank d — c and a matrix representation of full rank can be obtained by selecting a node in every 
connected component of G and deleting the corresponding rows. The matroid concepts above have 
natural interpretations in graph-theoretic terms: circuits correspond to cycles and independent sets 
to forests. The characteristic polynomial xdi) = XMci't) in this context is also called the tension 
polynomial and t^Xcit) counts the number of proper t-colorings of G where t € In the setting 
of Section [21 the relevant hyperplane arrangement is the graphic arrangement associated with G, 
which has the defining polynomial 

/g(z) = n ~ ^i)- 

The entropic discriminant i7G(b) is the equation of the branch locus of the gradient map V/^,. As 
Ag does not have full rank, we assume = for the rows i that were deleted when passing from 
Ag to a rank d — c matrix with d — c rows. The gradient map V/^ is discussed in [111 Remark 8]. 

Example 19 (Cycles). Let G = Gd+i be the cycle with n = d+\ edges. Every collection of d or 
fewer edges is independent and Mg has a unique circuit. The truncated matrix Aq+i has corank 
1 and Mg is the uniform matroid Ud^d+i- The reciprocal plane is a hypersurface of degree d, 
and the entropic discriminant Hg^^^Oo) is the polynomial of degree d{d — 1) seen in Section [3l o 
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Example 20 (Complete graphs). As the name says, the complete graph G = K^+i has all possible 
edges on d + 1 nodes. The characteristic polynomial is the chromatic polynomial divided by t: 

XK.^t) = {t-l){t-2)---{t-d). 
The reciprocal plane C,^^^^ is a projective variety of degree (— l)'^X-ft:d+i (0) = ^' • We find that 

degF,,,,(b) = ^\d-l-\-\-\...-\).dl 

is the value of the matroid invariant Q for the incidence matrix A^^ of the complete graph Kn- 
For example, for d = 3 we get the complete graph on 4 nodes, with node-edge incidence matrix 



1 ■ 

1/ 

The reciprocal plane C~^^ is a surface of degree 6 in P^. Its homogeneous prime ideal is generated 
by four quadrics, one for each of the 3-cycles in K4^. The entropic discriminant Haj^^ defines a 
curve in the projective plane P^. That curve has degree 14 and it has precisely six real points, o 

The matroids associated with the retina equations ([6]) are different from the matroids Mq above. 
Their matroids correspond to all-negative graphs in Zaslavsky's theory of signed graphs [28j. Here, 
an all- negative graph — G is an ordinary graph with all edges marked by —1. The incidence matrix 
A-G of —G has entries in {0, 1} where a 1 signifies an incident node-edge pair. The corresponding 
matroid M{—G) = M{A-.g) is the unoriented cycle matroid. The matroid-theoretic notions for 
M{—G) translate to (signed) graph concepts but the transitions are more involved. For all-negative 
graphs, the circuits correspond to even primitive walks, that is, even cycles or pairs of odd cycles 
connected by a simple path (of length possibly 0); cf. [281 Cor. 7D.3(e)]. For the state of the 
art on algebraic properties of the circuits of A-g see the recent work of Tatakis and Thoma |24j . 
Evaluations of the characteristic polynomial have interpretations in terms of signed colorings |27j . 

For example, the all-negative complete graph —K4 on four nodes has the incidence matrix 

1 1 0\ 
110 
10 10 1' 
10 11/ 

Note that this matrix has rank 4. Its matroid has the characteristic polynomial 

X-K^t) = - + 15^2 _ 17^ + 7. 

The characteristic polynomials for the all-negative complete graphs on any number of nodes were 
computed by Zaslavsky [271 Eqn. 5.8]. We presented his formula in the introduction in ([8|). An 
equivalent formula in terms of generating functions due to Stanley [22\ Ex. 5.25] was shown in Q. 

For the matrix (j28p . the reciprocal variety ^Z^^ is defined by the three cubic equations 

X12X13X2A - a;i2Xi3X34 - X12X2AX34 + Xi3X24a;34 = 0, 

a;i3a;i4a;23 - a;i3a;i4a;24 - a;i3a;23a::24 + a;i4a;23a;24 = 0, 
a;i2a;i4X23 - xi2XiiXu - a;i2a;23a;34 + xiiX2zX34: = 0. 

The task in ([B]) is to solve these cubic equations together with linear equations A^k^ • x = b for the 
six unknowns X12, • • • ,X34. The number of solutions to this system is ^{M-k^) = 7, and all seven 
solutions are real when the hi are real. One of the solutions has only positive coordinates if and 
only if the column vector (61, 62, ^3, ^4) of parameters lies in the convex polyhedral cone spanned by 



A 



1 


1 


1 








1 








1 


1 





-1 





-1 











-1 





-1 



(28) 



A. 



K4 



(\ 
1 
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the columns of A^m- The entropic discriminant H^x^ibi, 62, &3, ^4) characterizes parameter values 
for which the number of solutions is less than 7. It is a surface in of degree 2(4 ■ 7 — 17) = 22. 
The Mobius invariant ^{M^Kj) and the degree of H^k^ for larger values of d are displayed in ([TO]) . 

We close this section with the remark that the study of characteristic polynomials of matroids 
is an active area of current research in combinatorics. The coefficients of x(i) have interpretations 
as face numbers of broken circuit complexes and form a log-concave sequence. This log-concavity 
was a longstanding conjecture recently resolved by Huh ^JLj for graphs and in its full generality by 
Huh-Katz [12]. Their methods of proof are based on the geometry of reciprocal planes, our topic in 
the next section. Specifically, a key player in [12j is the tropicalization of the graph of £ --^ C~^. 



Entropic discriminants arise as branch loci from projecting reciprocal planes. This was already 
hinted at in the proof of Proposition HI We shall make this precise in Section [H where it will be 
our main ingredient in the proof of Theorem O In this section we build up to this proof by deriving 
some results on reciprocal planes. We believe that these results are of interest in their own right. 

We fix a d X n-matrix A of rank d with no zero columns. Its rows span a (d — l)-dimensional 
subspace C in the projective space P"~^. Let T denote the dense torus in P"~^, i.e. the complement 
of the n coordinate hyperplanes {xi = 0}. The reciprocal plane is the Zariski closure of the 
coordinate- wise inverse of £nT, as in ([5]). It is an irreducible projective variety of dimension d—\. 
The inversion map from C to is birational and it is an isomorphism on £ fl T. 

Proudfoot and Speyer [19] showed that is stratified by the flats of the matroid M[A). 
Recall that J C [n] is a flat of M{A) if and only if rk(Aj) < rk(^j/) for ah J' D J. Here 
A J denotes the column-induced submatrix of A. For a flat J C the corresponding stratum 
n P"^ = {p G : supp(p) C J} is isomorphic to £7^, the reciprocal plane associated to the 
restriction Aj. We shall investigate these boundary strata and the singular locus Sing(£~^) of £~^. 

Let C[A) C denote the set of circuits of M{A). To each v G C{A) we associate a polynomial 



These circuit polynomials cut out the variety £ . In fact, Proudfoot and Speyer proved the much 
stronger result that {hy : v G C{A)} is a universal Grobner basis for the prime ideal of £^^. 

As the set of all circuits is typically rather large, one might be interested in a smaller set of 
polynomials to cut out £~^. The following characterizes subsets of the set of circuit polynomials 
that cut out £~^ set-theoretically. As we saw above, the boundary of £~^\T in £~^ is described 
by flats of M{A). Recall that J Q [n] is a flat if and only if \J^ f] supp(t')| 7^ 1 for every circuit 
V G C{A). We say that a non-flat J C [n] is exposed by a circuit v G C{A) if \J^ n supp(t;)| = 1. 

Proposition 21. Let B C C{A) be a subset of the set of circuits. The corresponding set of circuit 
polynomials {h^ : v ^ B} cuts out £~^ set-theoretically if and only if B exposes every non-flat. 

Proof. We first prove the only-if direction. Suppose that J is a non-flat that is not exposed by any 
f G S. Then, for each v £ B, either | J*^ n supp(f)| > 2, in which case is identically zero on P*^, 
or \J^ n supp(?;)| = 0, in which case f is a circuit of Aj and hy vanishes on jOj^. This shows that 
the subvariety of P"~i cut out by {h^ : v €^ B} contains Cj^ and hence is strictly larger than £~^. 

For the if-direction we assume that B exposes every non-flat. Let p be any zero of {h^ : v B}, 
and let J = supp(p). Suppose that J is a non-flat of AI(A). Then there exists v £ B that 
exposes J. This means that exactly one of the terms of is non-zero at p, and hence h^{p) 7^ 0. 
We conclude that J is a flat of M{A). Since is a boundary stratum of £~^, it is sufficient 



5. Geometry of Reciprocal Planes 



(29) 
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to prove that pj^ G rowspaiiAj. For this, we shah prove that the kernel of Aj is spanned by 
{vj : V € B and supp(f) C J}. Let Jq C J be a basis of M{Aj). If Jq is not a flat, then there is 
a circuit vi € B supported on J such that supp(vi)\Jo = {ji}- Set Ji = Jq U {ji} and repeat the 
procedure. This process terminates after k = \ J\ — rk(Aj) = dimkerAj many steps. The matrix 
of the resulting circuits vi, . . . ,Vk is lower-triangular and hence gives a basis for ker(Aj). □ 

This previous result highlights the connection of our study to tropical geometry. 

Remark 22. Combining Proposition[2T]with the results of [26], we infer that a collection of circuits 
cuts out the reciprocal plane set-theoretically if and only if it constitutes a tropical basis for the 
tropicalization of the linear space £.. Yuster and Yu [26', Sect. 2.2] showed that different inclusion- 
minimal tropical bases for C need not have the same cardinality. Specifically, the uniform matroid 
U2.5 has inclusion-minimal tropical bases of size 5 and 6. From this we can infer that Proposition [2T] 
holds only set-theoretically and not in the ideal-theoretic or scheme-theoretic sense. 

Example 23. If the matroid M{A) is uniform, then the prime ideal of the reciprocal plane 
is minimally generated by ("^ ) polynomials of degree d. This can be seen as follows. The initial 
ideal of C^^ with respect to the reverse lexicographic term order is generated by the square- free 
monomials representing broken circuits. These where 1 < ii < • • • < id < — 1, 

so their number is ("^^)- By \26\ Lemma 5], the basic circuits obtained by adding the last element 
n form an inclusion-minimal tropical basis for C~^. Hence, by Remark 1221 the corresponding 
minimally cut out C~^. It follows that they form a minimal generating set for the ideal of o 

We now come to the main result in this section, namely, the characterization of the tangent 
cone of the reciprocal plane at any point. For the sake of convenience, we here identify the 
(d— l)-dimensional projective variety with the corresponding d-dimensional afhne variety in C". 

The tangent cone TCpX of a variety X C C" at a point p is the closure of the scheme of 
secant lines of X through p. For a polynomial / € C[x], the initial form in_i(/) is the non-zero 
homogeneous component of / of minimal degree. The tangent cone TCpX is defined by the ideal 

(30) /(TCpX) = (in_i(/(x + p)) : feI{X)). 

The following result shows that the tangent cone of at any point is reduced and irreducible. 
Here we use JC-a/j to denote the {d — rank(J))-dimensional linear space in C"'/Cj ~ C"~I'^L 

Theorem 24. Let A € M'^^" be a matrix of full rank d and let Cr^ be its reciprocal plane in C". 
For any point p € CT^ with support J, the tangent cone is isomorphic to the direct product 

(31) TCp£-i ^ Cj xC-y^, 

where "=" denotes the equality of affine schemes after a linear transformation in C". 

Proof. We inspect the initial forms of the circuit polynomials that define C~^. Let v G C{A) be 
a circuit with support C = supp(w) and circuit polynomial /it,(x) as in (p9|) . First suppose that 
C J . We write v = v' -\- v" where supp(v') = C fl J and supp(f ") = C\J . Then v" is a circuit of 
the matroid M{A/ J) obtained from M{A) by contraction at J. The following identity holds: 

K{-K + p) = -K^i-K + p) + (x-Fp)'^""' • V(x). 

Every term of :xp\'^hyi{K + p) has degree at least |C\J| while 

in_l((x+p)^n-^/l,.(x)) = p^""' -Kni-K) 

has degree |C\J| — 1. This means that hv"{x.) is the initial form of /i^(x + p). As every circuit w 
of the contraction M{A/J) is the restriction v" of some circuit v of M{A), we conclude that the 
tangent cone ideal at p contains the prime ideal (/iu;(x) : w G C{A/J)) that defines J~-a/j- 
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Next suppose that C J. Then p is a regular point on the hypersurface {h^ = 0}, and the 
initial form in_i(/i^(x + p)) is the differential Dph^. The differential of at the point p is 




The second equality holds because p ^ lies in rowspan(yl j) n (C*) , and the third equality follows 
from the fact that —Vi/pi = Ylj^i'^i /Pj-> since t> is a circuit for Aj. Thus, Dph^ vanishes on the 
rowspan of Ajdiag(pj)^, denoted Cj{p), and all circuits vanishing on this row span arise this way. 
We have shown that the prime ideal of the irreducible variety Cj{p) x is contained in the 

ideal of the tangent cone of at p. Since both ideals have the same height, and the former is 
prime, it follows that they are equal. This proves the equality of schemes that was claimed. □ 

A closer inspection of the proof reveals that the initial forms of /i^(x + p) for v G C{A) furnish 
a universal Grobner basis for the tangent cone of at p. In particular, we obtain a simple 
description of the tangent space of at a point p by taking those initial forms that are linear. 

Corollary 25. For a point p E C^^ with support J, the tangent space is orthogonal to the space 
spanned by the circuits of the d x \ J\-matrix Aj diag{pj)'^ and the circuits of A/ J of size 2. 

Proof. The tangent space is cut out by the linear forms in the ideal of the tangent cone. From the 
initial forms in the proof of Theorem 1241 we see that in_i(/i^) is linear whenever |supp(u) nJ'^\ < 2. 
If C C J, then in_i(/it,) corresponds to a circuit of Aj diag{pj)^ . Otherwise, the two elements of 
C\J are parallel in the contraction A/, J and the corresponding circuit polynomial is linear. □ 

Our discussion shows that the dimension of the tangent space is constant on each stratum of 
C^^. We obtain the following characterization of the singular locus of the reciprocal plane C~^. 

Corollary 26. The singular locus of the reciprocal plane C^^ is pure of codimension 2. It is the 
union of all boundary strata Cj^ such that the contraction M{A/J) is a non-basic matroid. 

Proof. A point p E C^^ is smooth if and only if the codimension of the tangent space equals 
codim(>C~^) = n — d. The description of the tangent space in terms of the matroids M{Aj) and 
M{A/J) in Corollary[25]shows that its codimension is | J|— rk^j+Par(j4/J) where Par(A/J) is the 
dimension of the space of 2-circuits of A/ J. Suppose M{A/J) has r distinct 1-flats (or lines), and let 
Ai , . . . , be the sizes of these parallelism classes. The circuits of each parallelism class span a linear 
space of dimension Aj — 1. As these circuits are disjoint, we have Par(A/J) = X]I=i(-^«~l) — \ J'^\—i'- 
The number of parallelism classes of M(A/J) is at least rk(^/J). Thus the codimension of the 
tangent space is < n — d, and equality holds if and only if M{A/J) is basic (cf. Proposition fTT]) . 

Finally, to see that the singular locus is pure of codimension 2, we note that if M is any non-basic 
matroid if rank r > 3, then there is an element e such that M/e is non-basic. To show this, we can 
assume that M is non-basic on r + 1 elements, each representing a different line. If M = Mi ©M2 is 
not connected and Mi is non-basic then any e G M2 will work. Otherwise, M is a uniform matroid 
and the contraction is uniform of rank r — 1 > 2 on r elements. By Example [TSl the uniform 
matroid Un^d is non-basic if and only if n > d > 1. Therefore, if J is a flat of M(A) such that 
M{A/J) is non-basic of rank > 3, then there is a flat J' D J such that M{A/J') is non-basic. □ 

6. Ramification Locus 

The entropic discriminant describes the locus of points b € W^^^ such that the zero-dimensional 
scheme defined by the constraints x G and Ax = b is not reduced. Equivalently, the entropic 



16 



RAMAN SANYAL, BERND STURMFELS, AND CYNTHIA VINZANT 



discriminant is the defining polynomial of the branch locus of the map A: C"^ ^ P'^-^ We be gm 
with the observation that this map has no base points and is hence a projective morphism. 

Lemma 27. The variety is disjoint from the center of the projection A : P"'"^. 

Proof. Our claim states that r\kev{A) = {0} holds in C". Let phe a vector in C^^ r\ker{A) and 
J = supp(p). Then pj^ = z Aj for some z € C"', and p G ker(A) implies = z{Ap) = {zA)p = 
J2jeJ Pj~^Pj ~ l"^!- follows that J = % and p = 0. □ 

We now focus on the ramification locus of the dominant projective morphism A : C^^ — >■ P'^^-'^. 
By definition, this is the Zariski closure of the set of regular points p € for which 

(32) C + rowspan Jac(/:~^)(p) / C". 

Here Jac(jC~^) is the Jacobian matrix of C~^, whose row vectors are V/it,(x) for v G C{A) in (|29p . 
This condition states that the intersection of and {x : ^x = Ap} is not transverse at p. 

The ramification scheme TZa = Proj(C[x]/ Ja) is defined by the following ideal in C[xi, . . . , 

(33) Ja = + (ny.n minors of ^j^^^^-ij^ )) • ( x (n—d) minors of Jac(£^"'^) 

By the Zariski-Nagata Purity Theorem [16J, the ramification locus is pure of codimension 1 in 
C^^. Hence the ramification scheme TZa is either empty or has codimension 1 in C~^. The former 
happens when A is basic, and the latter happens when A is non-basic. We prove in Section [7] that 
TZa contains the singular locus of and hence that the saturation step in (p3]) is redundant. 

Definition 28. Let A G R*^^" be a non-basic matrix of rank d. The ramification cycle is the 
algebraic cycle of dimension d — 2 in p*^"^ defined by the ramification scheme TZa- By Corollary [U 
the push-forward of the ramification cycle under the morphism A : — >• P'^"^ is a cycle of 
codimension 1. We define the entropic discriminant of A to be the homogeneous polynomial HaO^) 
that represents this cycle in P'^"^. It is unique up to multiplication by a non-zero constant. 

The following example shows that the ramification cycle may not be reduced. 
Example 29. Let A be the matrix in Example El For a ^ 0,2,3, the prime ideal of C^^ equals 
I{C^^) = { 2x\X2 — 3x13:3 -|- X2X3, 2x\X2 — ax\X!^ -|- (a — 2)x2a;4, 3x1X3 — 0x1X4 -|- (a — 3)x3X4 ). 
The ramification ideal J a is the sum of I{C^^) and the ideal of 4 x 4 minors of the matrix 



A 

jac(/:-i; 



/II 1 1 \ 

2 3 a 

2x2 — 2xi + X3 — 3xi + X2 

2x2 ~ axi 2xi — (a — 2)x4 — axi + (a — 2)x2 

8x3 — ax4 3xi + (a — 3)x4 — axi + (a — 3)x3 

y X3 — (a — 2)x4 X2 + (a — 3)x4 —(a — 2)x2 + (a — 3)x3 J 



The ramification cycle is a zero-dimensional cycle of degree 4 in P^. For the special value a = 6, it 
is twice the reduced cycle of degree 2 defined by (2x2 — 3x3 + 6x4, 2xi — X3 -|- 4x4, x\ — 4x3X4 -|- 8x|). 
The push- forward of this cycle under P^ P^ is defined by the binary quartic in Example [3l o 



Since the projection A : P""i --^ pd-i j-^g^g base points on the subscheme TZa (by Lemma [2l 
the push-forward by A preserves the degree of the ramification cycle. Thus, in order to establish the 
degree formula in Theorem^ it suffices to show that the degree oiTZA equals 2(— l)'^((ix(0)-|-x'(0)). 
In order to compute its degree, we use a slightly different description of TZa- Let T denote the 
dense torus {xiX2 • • • x„ / 0} in the projective space P"~^. Inside T, the variety is a complete 
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intersection. Namely, it is defined by -B • x ^ = 0, where B = {Bi, . . . , is a Gale transform for 
A, that is, an [n — d) x n-matrix whose rows span the kernel of A. Consider the polynomial 

i=l 

The n x n-matrix above now plays the same role as the Jacobian matrix did in (j33p . Thus the 
hypersurface defined by ^^(x) = inside fl T is the restricted ramification locus IZa H T. 

If QA is zero at a point p £ T then the intersection ker(j4 diag(p)^) fl ker(i?) contains a non-zero 
vector. The kernel of B is spanned by the rows of A, so the d x d-matrix ^diag(p)^yl"^ also drops 
rank. Hence 5a(x) divides det(74 diag(x)^A-^). Both polynomials have the same degree 2d, and 
hence they are equal (up to a scalar). Using the Cauchy-Binet Formula, this gives 

(34) 5a(x) = det(^diag(x)2^^) = ^ det{Ajfl[xl 

We next define similar polynomials that cut out TZa on the non-singular boundary strata of 
Let J C [n] be any proper fiat of rank r in M{A) and set xj = (xj : j (z J). Let Aj now denote 
any r x | J| submatrix of Aj = {Aj : j (z J) whose rows are linearly independent. We define 

(35) 9AjM = det(ijdiag(xj)2ij) = ^ det(i/)' J] xf . 

Hi) 

Here Ai denotes the square submatrix of Aj induced on the r columns indexed by / C J. 

Lemma 30. Let p be a smooth point on the reciprocal plane Cr^ with supp(p) = J. Then the 
ramification locus TZa contains the point p if and only if gAj{pj) = 0. 

Proof. Since p is smooth, the condition ([32]) reduces to ker{A j) n ker(Jac(/3j^)) ^ {0}. From the 
argument prior to we see that, for G (C*)'^, this is equivalent to gAj{pj) = 0. □ 

Remark 31. This characterization shows that the ramification locus IZa equals the closure of 
its intersection with the torus, TZa H T. To see this, suppose that TZa has some component Z 
contained in the boundary of the torus {xi • • • x„ = 0}. Then Z is contained in for some 
proper fiat J, where dim(>Cj^) = rank(ylj) — 1. Since TZa is pure of codimension-1 in C~^, we see 
that dim(Z) = d — 2. It follows that rank(Aj) = d — 1 and Z = jOj^. However, M{A/ J) has rank 
1 and is therefore basic. Lemma [30l then tells us that is not contained in TZa- This shows that 
to define the ideal Ja in (p3]) . we could instead saturate with respect to the ideal (xiX2 • • • x^). 

We shall now use the polynomial 5^(x) to compute the degree of the ramification cycle. 

Proof of Theorem\^ Let A be a non-basic real d x n-matrix of rank d and xi^) the characteristic 
polynomial of the matroid M{A). We shall prove that the degree of the algebraic cycle underlying 
the {d — 2)-dimensional subscheme TZa of P"~^ equals the matroid invariant Lemma [271 then 
implies that Ha{^) has the same degree, and this will complete the proof of Theorem [2j 

From above, we know that the scheme TZa is contained in the hypersurface {qa = 0} of P"~^. 
Let TZa denote the scheme-theoretic intersection of the reciprocal plane with this hypersurface: 

(36) Ua = Proj(C[x]/(/(/:-i) + (5A))). 

The (d — 2)-dimensional scheme TZa is the intersection of the {d — l)-dimensional irreducible variety 
and the hypersurface qa- By Bezout's Theorem [9; Thm. 1.4.4], its degree equals 

(37) deg(^^) = deg(5A)-deg(/:-i). 
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We claim that TZa decomposes into #Hyp(A) + 1 components of dimension d — 2, one of which is 
TZa- Here Hyp(A) denotes the set of hyperplane flats, that is, flats J such that rk(^j) = d—1. We 
see that TZa and TZa agree in the torus T, so their difference must lie in the coordinate hyperplanes. 
Recall from Section [5] that the reciprocal plane intersects the dense torus T"^ of P"^ if and only if J 
is a flat, and if so, the closure of that intersection is the reciprocal plane jCj^- Such a stratum has 
dimension d — 2 in P"~i if and only if J is a hyperplane flat. Since J G Hyp(A) does not contain a 
basis of M{A), each summand in the formula (I34p for gA vanishes on T"^ . To be precise, gA vanishes 
to order exactly 2 on the torus T'^ , since J is only one element away from containing a basis. 

Furthermore, the strata Cj^ are not contained in TZa for J S Hyp(A). This follows from 
Lemma [30l Indeed, by Corollarv 1261 the points in n T'^ are non-singular in C~^, and hence 
the polynomial gAj{^j) is not identically zero on Cj^. We conclude that the irreducible varieties 
jCj^, for J G Hyp (A), are components of dimension d — 2 and multiplicity 2 in the scheme TZa- 

We have derived the following equidimensional decomposition of the cycle defined in (|36p: 

(38) TZa = TZa U I \J 2 ■ Cf 

\jeHyp(A) 

Since the degree is additive on equidimensional cycles, we can use ([37|) to conclude that 

(39) deg(7^A) = deg(9A)-deg(/:-i) - 2 ^ deg{Cf) = 2d • /x(A) - 2 ^ ^(^j). 

JgHyp(A) JeHyp(A) 

The coefficient oft* in the characteristic polynomial x(t) equals (— 1)*^"* times the sum of the Mobius 
invariants fJ.{Aj) where J runs over all flats of rank d — i. For i = this gives ij.{A) = (— l)'^x(0)i 
and for i = 1 we get Ylje]iyp{A)l-^i^j) ~ (~l)'^^^x'(0). Hence the right hand side of (p9]) equals 
the desired matroid invariant (j4]) . This completes the proof of Theorem [2j □ 



The decomposition (|38p can be used to compute the ideal of the ramification scheme. Namely, 
since all hyperplane strata Cj^ lie in complement of the torus T, we have the algebraic identity 

(40) Ja = {l{C'') + {gA)):{xiX2---Xnr. 

We illustrate the identity (|40p and our proof of Theorem [2] for the codimension 1 case. 

Example 32. Let A be the matrix in equation (|2Up of Section [31 The reciprocal plane is the 
hypersurface defined by the elementary symmetric polynomial e„_i(xi,X2, . . . ,Xn). The equation 
([3ip defining the ramification locus in the torus is gA = en-i{x\,x\, . . . ,x^). The scheme TZa in 
()36p is the complete intersection of these two hyper surf aces. Its ideal has the primary decomposition 

, , ^en_i(xi, X2, • • • , Xn)^ (in—l{Xi, 2^2' • • • ; "^n)) 

= (en-i(2;i , X2 , • • • , Xn ) , Gn—2 

(Xl , X2 , • • • , Xn 

)) n n 

This is the decomposition discussed after (j37p . with the first intersectand being the ideal J a that 
defines TZa- This ideal is contained in the Jacobian ideal of the reciprocal plane because 

1 ^ den-i 

2 = 1 

This identity proves the ideal-theoretical inclusion J a C I(Sing(i2~^)). We conclude that TZa 
contains Sing(£~^) when n = d + I. As we shall see in Section [3 the inclusion Sing(£~^) C TZa 
implies that the real variety of HaO^) is indeed the union of codimension 2 planes given in ()27p . o 



We close this section with a combinatorial proof of the assertion, stated informally immediately 
after Theorem [2l that generic matrices maximize the degree of the entropic discriminant. 
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Proposition 33. ^45 A ranges over all non-basic d x n-matrices of rank d, the degree of the entropic 
discriminant ff^(b) attains its maximal value 2{n — d)(^Z2) when the matroid M{A) is uniform. 

Proof. It follows from Theorem [2] and Example 1181 that 2(n — c^)(]^Z2) degree of the entropic 

discriminant when M = M{A) is uniform. We must show that this number is a strict upper bound 
otherwise. The claim is an entirely matroid-theoretic statement, and so let us define 5{M) = 
2(— l)'"'^(*^)(rk(M)xM(0) + Xm(^)) matroids M. The characteristic polynomial satisfies a 

deletion-contraction recurrence, namely, XM{t) = XM\ei't) ~ XM/ei^) for e S M not an isthmus. It 
follows that the entropic degree satisfies a deletion-contraction recurrence plus a correction term: 

5{M) = 5{M\e) + 5{M/e) + /u(M/e). 

All three terms on the right hand side are non-negative. The desired inequality follows by induction 
on the rank d and corank n — d. In rank 1 all simple matroids are uniform. Corank 1 is dealt with in 
Section[3j The same argument shows that n{M) < (^Zi) with equality if and only if M = Ud^n- D 

7. Real Issues 

Our point of departure for this paper was the observation that, for real b, the equations ([2]) have 
only real solutions, namely, the ii{A) analytic centers of the bounded regions in the arrangement 
of n coordinate hyperplanes in {Ax = b} ~ M""*^. It is thus natural to ask what it means for two 
such analytic centers to collide, and how this relates to the real points in the ramification locus 
and in the entropic discriminant. We shall prove that the real loci of these two complex varieties 
are both pure of codimension one. Our first step in this direction is the following lemma. 

Lemma 34. All real points in the ramification scheme are singular in the reciprocal plane: 

(42) (7^A)R C Smg{C-%. 

Proof. The sum of squares formula in ()34p reveals that ^^(x) = has no real solutions in the torus 
T. In symbols, (7^^nr)]R = 0. Likewise, for any flat J with Cj^ nonsingular in £~^, the polynomial 
gAj is a similar sum of squares, and hence (T^yinT'^)^ = 0. Lemma [30] ensures that no regular point 
of with real coordinates lies in the ramification locus of the morphism A : Cr^ — ?• P'^^^. □ 

The following is our main result in this section. We find that the reverse inclusion holds in (|42p . 
Theorem 35. The ramification scheme TZa contains the singular locus of , and we have 

(43) (7^A)R = Sing(/:-i)M. 

This theorem implies that the saturation in the formula (j33p for the ramification ideal J a was 
unnecessary. Before presenting the proof, we shall derive two corollaries and discuss one example. 

Corollary 36. The Zariski closure of {JZa)r is pure of codimension 2 in . 

Proof. Theorem [35] implies that the Zariski closure of the set (7^a)r of real ramification points 
equals the singular locus Sing(£~^) of the reciprocal plane . Corollary [26] represents Sing(>C~^) 
as a union of linear spaces all of which are defined over M and have codimension 2 in . □ 

We now obtain the following characterization of the real locus of the entropic discriminant. 

Corollary 37. The Zariski closure of the set of real points in the hypersurface defined by the 
entropic discriminant HaO^) is pure of codimension 2 in P"'^"'^. Its irreducible components are the 
linear spaces span{Aj : j G J), where J runs over all non-basic corank 2 fiats of M{A). 
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Proof. The real variety of is the image of the real points in TZa under the fi{A)-to-one morphism 
A : — > P*^"^. Hence the real variety of Ha is pure of codimension 2 in F'^-^ as well. The 
description of its irreducible components now follows from that given in Corollary [26l □ 

We now revisit our very first example to illustrate the previous corollary. 

Example 38. For d = 3, the codimension-2 strata of are the n coordinate points Cj in P"~^. 
Their images under the map A are the columns Ai, . . . , An- For generic A, the points ei, . . . , e„ 
comprise Sing(£^^). Lemma [34l implies that Vs_{Ha) is contained in {Ai, . . . , An}, and Theorem [351 
reveals that equality holds. For special 3 x n-matrices A, the matroid M{A/i) may be basic for 
some i. If this happens then Cj is a non-singular point in and its image Ai does not belong to 
Vr{Ha)- Looking back at Exampledl we notice that the matroid M{A/i) is basic for i = 1 and it 
is non-basic for i = 2,3, 4, 5. This explains our finding in ([3]) that the real variety Vm.{Ha) consists 
of precisely the four points A2, A3, and A^ in the projective plane P^. o 

We are now ready to present the proof of our main result in this section. 

Proof of Theorem \35[ We first note that the identity (j43p follows immediately from Lemma [M] and 
the inclusion TZa 5 Sing(>C~^) in the first assertion. Hence it suffices to prove that inclusion. 

By Corollary [26} the singular locus of is a reducible variety whose irreducible components 
are the boundary strata Cj^ where M(A/J) is a non-basic matroid of rank 2. We consider one 
such component Cj^, regarded subvariety of C'' x {0} inside of C" = C"^ x C"''. A generic 
point of Cj^ has the form {p,0) where p £ (C*)"^. Our goal is to show that this point lies in the 
ramification locus TZa by producing a sequence of points in TZa that converges to (p, 0). 

We may assume that J = {1, . . . , A;} is a fiat of rank d — 2 and our matrix A has the block form 



A 


* 





B 



where A G M('^-2)x'^' and B G are both of full row-rank. In these coordinates, we get 

M{Aj) = M{A) and M{A/J) = M{B). 

Now, let us return to our generic point (p, 0) G jC-j^- The partial specialization (^^(p, xjc) is a 
polynomial in C[xk+i,Xk+2, ■ ■ ■ ,Xn]- It is non-homogeneous and its terms of lowest total degree 
come from those bases / of M(A) for which |/ H J| = d — 2. For any such I, we have 

det{Ai) = det(i/nj) ■ det{Bjnjc). 
From this we see that the initial form of gA{p,^j'') of lowest degree terms can be written as 

(44) in_i (c/a(p,xjc)) = • 

From the results of Section[6]we know that {g^ = 0}n Cj^ has codimension 1 in Cj^- This implies 
g^{p) 7^ because the point (p, 0) was chosen to be generic in Cj^. 

In order to proceed, we need to represent the ramification locus around p by a single polynomial, 
rather than as a subvariety of C~^. To do this, we rationally parametrize the points xjc for which 
(p, xjc) lies in using the matrix B. First, note that the intersection of the linear space C with 
{p} X C'^'' gives an affine linear space in C'^'' of the form v + rowspan(i?) for some vector v in C'^''. 
We can parametrize this space by v + zB where z = (zi, 2:2). This gives the rational parametrization 
(p, {v + zB)~^) of the intersection of with {p} x C'^" . 

Now we plug this parametrization into ^(^(p, xjc) and clear denominators to get a polynomial in 
C[zi, Z2]. Define g{z) G C[zi, Z2] to be this polynomial, 

(45) g{z) = gA{p,{v + zB)-')l[{v, + zBif = ^ det{Aif J] pj J] (vj+zBj)^. 
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If z is a solution to g{z) = for which each coordinate of v + zB is non-zero, then the point 
{p, {v + zB)~^) lies in the ramification locus TZa- 

Since J is a flat, the n — k linear forms zBi are non-zero for all indices i. This implies that 
Xi = + zBi) has degree —1. Thus the terms of highest degree in g{z) correspond exactly to 
the terms of lowest degree in (//^(p, xjc). From (I44p . we see that the leading form of g{z) is 

ini(5(z)) = gj^{p)-gB{{zB)-^)-\{{zB,f. 

Our next step is to find a solution to the initial equation inig'(z) = and to then extend it to the 
desired sequence of points in TZa- As the matroid M{A/J) = M{B) is non-basic, it follows from 
Corollary [5] that the ramification TZa/j is nonempty. Hence there is a point q G Ca/j H (C*)"~*^ 
such that gA/j{Q^^) = 9b{q^^) = 0. Let z be the unique vector such that zB = q. We may assume 
that B has the form (Id2 B' ). Thus implying that Zi = qi ^ for i = 1,2. 

By Lemma HOl below, we can extend this solution z € (C*)^ to a solution Z = Z(e) of g{z), where 
the coordinates of Z = (Zi, Z2) lie in the field C{{e}} of Puiseux series: 

Zi = Zi- + higher order terms S C{{e}} fori =1,2. 

Moreover, by Lemmas 1391 and 1401 these series converge in a neighborhood of zero in M>o. 

Now consider the point Q = Q{e) = v + ZB with coordinates Qi = qi^ + ■ ■ ■ in C{{e}}. We can 
invert Qi in the field of Puiseux series to get 

= Q^^e + higher order terms G C{{e}} for i = 1, . . . , n — /c, 

and these series converge for real e in an open segment (0,eo) near zero (see Lemma [5^ below) . 

Then, by (05]), the point {p, Q-^) in C'^ (S)cC{{e}} is a zero of the polynomial ^'^^(x). Specializing 
to sufficiently small e G M>o, gives a point {p,Q{e)~^) € (C*)" that belongs to the ramification 
locus TZa- Furthermore, as e approaches 0, the limit of the points {p, Q~^{e)) is {p, 0) in Cj^ x {0}. 
This shows Cj^ C TZa and consequently Sing(£~^) C TZa- D 

Before Lemma HOj we need a short lemma on the convergence of reciprocals of Puiseux series. 

Lemma 39. If x{e) is a nonzero Puiseux series that converges for e > in a neighborhood ofO, 
then its inverse x{e)~'^ in C{{e}} also converges for real e in an open segment (0, eo). 

Proof. Suppose x{e) = ue^ + higher order terms. We can write the field of Puiseux series as the 
union of C((e^/™')) over m G Z+. Thus for some m € replacing e with yields a Laurent 
series x(e™'), which also converges in a neighborhood of 0. In particular, e~^^x{e"^) is a convergent 
power series with constant term u and has an inverse y{e) in the ring of convergent power series 

(see P §6.4]). Then y(e) = l/u-\ satisfies e~™'^x(e'")y(e) = 1. Replacing e with e^/™, we see that 

€~^y{e^^™') is an inverse for x(e). Furthermore, since y(e) and y{e^^"^) converge in a neighborhood 
of zero, x(e)^^ = e~^y{e^/'^) also converges for e > in a neighborhood of zero. □ 

Now all that remains is to lift roots of initial forms to solutions over C{{e}}. 

Lemma 40. Let g{zi,Z2) be a polynomial with complex coefficients and initial form mi{g), con- 
sisting of the highest terms with respect to total degree. Let u = (mi,M2) € (C*)^ be any solution 
to the equation ini((^)(ui, U2) = 0. Then there exists a vector v{e) that satisfies g{v{e)) = and 
whose coordinates are Puiseux series of the form 

Vi{e) = Ui— + higher order terms in e, for i = 1,2, 

that converge for e in some neighborhood (0, eo) of zero. 
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Proof. We invert the variables Zi and work with the polynomial 

g{z) = .^(^)..^(^).5(^rS^2"')- 

The highest-degree terms of g then correspond to the lowest-degree terms of 'g. Furthermore, the 
point = {1/ui, I/U2) is a solution of in_i(5). 

Our hypothesis states that the Newton polygon of 5(2:1, 2:2) has an edge of slope —1, and 
(l/ui, l/ti2) is a root of the corresponding binary form in_i(5)(zi, Z2). Using the classical Newton- 
Puiseux algorithm, we can construct a power series expansion of Z2 in terms of zi = -^e, having 
the form Z2 = + ■ ■ ■ ■ The resulting series in e converges by the arguments in [8j §7.11]. 

This solution has an inverse in the field of Puiseux series, and this inverse will be our desired 
solution (fi(e), W2(e)) of 5(^1,22) = 0. Namely, if w{€) = {wi{(),W2{e)) € C{{e}}^ is the solution to 
^(zi, Z2) found in the paragraph above, then '^^(e) = Wi{e)~'^ is a solution to g{zi., ^2)- By Lemma [39} 
the Puiseux series Vi{e) converge in a neighborhood (0, eo) of the origin in M>o. □ 

This concludes our study of the entropic discriminant. The discussion above shows that there 
are still many unresolved problems concerning ff^(b). We list four open questions: 

Open Questions: 

(1) Is the entropic discriminant -ffyi(b) always a sum of squares? 

We know that the answer is yes for n = d + 1 and for d = 2, but even the case d = 3 of 
plane curves is open. It would be especially nice to write H^ih) as sum of squares in the 
maximal minors of the matrix {A,h), as we did in (jl6p and (jl7p for {d,n) = (2,3), (2,4). 

(2) What is the Newton polytope of the entropic discriminant H^ih)? 

For instance, when A is the matrix in (j20p then the table (|26p suggests that the Newton 
polytope of -^^(b) is the standard permutohedron, scaled by a factor of two. 

(3) Find ideal generators for the ramification scheme. 

Here is a concrete conjecture about minimal generators of the ideal J a in (|40p . Fix n > d+2 
and a dx n-matrix A whose matroid is uniform. We know from Example 1231 that I{C~'^) is 
minimally generated by (^~^ ) polynomials of degree d. We conjecture that Ja has precisely 
C^^^) additional minimal generators of degree 2d — 2, namely, the restrictions to of the 
rational functions gA{^) / XiXj for some i,j G [n]. We can show that these rational functions 
are polynomials on and that they vanish on Sing(/3~^). Do they generate our ideal? 

(4) How is the entropic discriminant related to the Gauss curve of the central curve ? 

The degree formula for the Gauss curve in [3] §5] is essentially the same as the degree 
formula we derived for HA(i>)- What is the most natural geometric explanation for this? 
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